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Abstract. We prove qualitatively sharp heat kernel bounds in the setting of Fourier-Bessel ex- 
pansions when the associated type parameter v is half-integer. Moreover, still for half-integer v, 
we also obtain sharp estimates of all kernels subordinated to the heat kernel. Analogous es- 
timates for general v > — 1 are conjectured. Some consequences concerning the related heat 
semigroup maximal operator are discussed. 



1. Introduction 

Let J u denote the Bessel function of the first kind and order v and let \\ n ,v : n > 1} be the 
sequence of successive positive zeros of J v in increasing order. It is well known that for each 
fixed v > — 1 the functions 

4>n( x ) = d n,u X~ U ~ 1/2 (\ ntU x) 1/2 J u (\n,vX), Tl = 1, 2, ... , 

form an orthonormal basis in L 2 ((0, 1), dfi v ); here d UjV = \/2| X n (u Jiy+i(A n ^)| _1 are normalizing 
constants and the measure is given by 

dnv(x) = x 2v+l dx. 

The system {0^ : n = 1,2,...} is usually referred to as the Fourier-Bessel system. Another 
Fourier-Bessel system arises by considering the functions 

<(x)=^ +1 / 2 <(x), n=l,2,..., 

which form an orthonormal basis in L 2 ((0, 1), dx). Fourier-Bessel expansions are of interest 
and exist in the literature for a long time, see [281 Chapter XVIII]. The related convergence 
problems were investigated in [21 El [29] . For the study of several fundamental harmonic analysis 
operators in the Fourier-Bessel context see [I0l CEH H21 [HI [HI [HI [22] and references therein. 
Some interesting open questions concerning harmonic analysis of Fourier-Bessel expansions can 
be found in a recent paper by Betancor [JJ. 

Given v > — 1, consider the differential operator 

L„ = —A — — 1 d 



eigenfunction of L v with the corresponding eigenvalue A 2 h 

v2 



x dx ' 

which is symmetric and nonnegative on C 2 (0, 1) C L 2 ((0, 1), d/j, u ). Each 0^, n = 1, 2, . . ., is an 

L v 4>n = ^n, 

This leads to a natural self-adjoint extension of L u given by 

oo 
n=l 
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on the domain Dom£j, consisting of all / G L 2 ((0, 1), <iu„) for which the series converges in 
L 2 ((0, 1), cZ//„); here (f,4>n)d^ v is the nth Fourier-Bessel coefficient of /. Clearly, the spectral 
decomposition of L v is given by the <\> v n . 

The aim of this paper is to find sharp estimates for the kernels 

oo 

(1) Gr(x,y) = ^exp(-U^)<(x)C(y), x,y € (0,1), t > 0, 

n=i 

where a £ (0, 2] is the index of subordination. These are precisely the integral kernels of the 
semigroups 

W' a }t>o :={exp(-tC /2 )} t>0 

generated either by C u itself (a = 2) or by its fractional powers (a < 2). Needless to say, the 
cases a = 2, 1 are of prior importance since then ([1]) becomes the heat or the Poisson kernel, 
respectively, related to the system {4>n}- The behavior of G^' a (x,y) does not seem to have 
been studied before. Note that for short times t a direct analytic treatment of the series in 
(PQ) is a complicated matter except for a few particular cases discussed in Section [H First of 
all, the series defining G"' a (x,y) is highly oscillating, and its behavior is hidden behind subtle 
cancellations between the oscillations. Moreover, J u is in general a transcendental function and 
there is no explicit formula for its zeros. Finally, the power a of the eigenvalues in the argument 
of the exponential makes the situation even more sophisticated. These obstacles can be better 
understood by means of the explicit computations in Section [J] for the cases v = ±1/2, a = 1, 2. 

Our method of estimating G^' a (x,y) is based on a connection between the Fourier-Bessel 
context and the classical setting related to multi-dimensional Euclidean balls. More precisely, we 
show that the fundamental solution Q^' 2 (x, y) of the classical heat equation in the d-dimensional 
unit ball, subject to the Dirichlet boundary condition, is related to G v t ' (x,y) with v = d/2 — 1. 
This allows us to transfer known bounds for C/^' 2 (x,y) to the Fourier-Bessel framework on the 
interval (0, 1) and conclude qualitatively sharp estimates for G^ ,2 (x, y) when v is half-integer, see 
Theorem l3.31 Essentially the same procedure applies to the subordinated kernels Q^' a (x, y) and 

G v t ,a (x, y), see Theorem 13.41 We emphasize that the sharp bounds for Q^' a (x, y), < a < 2, are 
quite strong results and seem to have been obtained in a complete form only recently \30\ 126] ; 
this may be a bit surprising, in view of simplicity of the geometry of the underlying domain, 
the Euclidean ball. For large values of t sharp estimates of G^' a (x,y) can be derived directly 
from the series representation ([I]), thanks to the exponential decay. In this case all v > — 1 are 
covered, see Theorem 13.71 We remark that all the results can be immediately translated to the 
setting of the system {ip^} since the relevant kernels coincide up to the factor (xy) u+1 ^ 2 . 

The paper is organized as follows. In Section [2] we first discuss the connection between the 
Fourier-Bessel setting in the interval (0,1) and the situation associated with Euclidean balls, and 

then relate the kernels G^' a (x,y) and q[ ,a (x, y). In Section [3] we invoke necessary bounds for 

5 t (x, y) and transfer them to the Fourier-Bessel setting in the interval. We also establish the 
long time behavior of G v t ,a {x, y) by analyzing the defining series. Finally, Section [5] is devoted to 
various comments and remarks. These concern, in particular, mapping properties of the maximal 
operators of the semigroups {T^ ,a } t >Q-, < a < 2. 

Throughout the paper we use a standard notation. While writing estimates, we will use the 
notation X < Y to indicate that X < CY with a positive constant C independent of significant 
quantities. We shall write X ~ Y when simultaneously X < Y and Y < X. 

Acknowledgment. This research was started during the sojourn in Wroclaw of the second- 
named author during the winter of 2011. She acknowledges a financial support received from 
Politechnika Wroclawska and wants to thank for the hospitality offered. 
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2. Connection with multi-dimensional balls 

For d > 1, let B d = {x G R d : |x| < 1} be the unit ball in R d and denote S ,d_1 = <9£ d . It 
is well known (see for instance |19t Chapter 2, H]) that there exists an orthonormal basis of 
eigenfunctions associated with the Dirichlet Laplacian in B d . These are expressed explicitly by 
the functions (j)^ and spherical harmonics, as described below. For more details on spherical 
harmonics and their connections with symmetry properties of the Fourier transform we refer to 
[271 Chapter IV]. 

Given k > 0, let Vk be the space of homogeneous polynomials of degree k in W 1 . Taking 
harmonic polynomials in Vk and restricting them to S d ~ l we obtain the space 

H k = {P\ S d-i :PeV k and AP = 0}. 

The elements of Hk are called spherical harmonics of degree k. The dimension of Hk is finite, in 
fact we have 

fc = dimfffc = (2k + d - 2) ( ^,|/_~^ ! , d>2. 

The case d = 1 is degenerated: S° consists of two points and so Do = f i = 1 and dk = for 
k > 1. 

Let {Y^ : 1 < m < t)fc} be an orthonormal basis for Hk in L 2 (S d ~ l , da) , where a is the 
standard (non-normalized) surface area measure on Define the functions 

<L«=^ +d/2 ~ 1 (W)^(^|), n>l, fc>0, l<m<c) fc ; 
here and later on we tacitly assume that for d = 1 only k = 0, 1 are considered. The system 

{*$, m :n>l,*>0,l<m<?> fc } 

is orthonormal and complete in L 2 (B d , dx). Moreover, it consists of eigenfunctions of the Dirich- 
let Laplacian in B d , 



n,k,m n,k+d/2—l ri,fc,m' n,k,m 



= 0. 



Thus the associated heat kernel and the subordinated kernels are expressed as 

n>l fc>0 l<m<0fe 

where a £ (0, 2] is the subordination index. These are the integral kernels of the semigroups 

{r t (d) n t> o -{e* P (-*(-A)«/ 2 )h >0 , 
where —A is understood to be the nonnegative self-adjoint operator in L 2 (B d , dx) whose spectral 
resolution is given by the . The domain of this operator can be identified with the Sobolev 

space H^(B d ), cf. [I]. 

We now observe that the analysis of the radial case in the context of expansions with respect 
to m } reduces to the Fourier-Bessel setting in the interval (0, 1) with the type index v = dj 
2 — 1. Indeed, if F(x) = /(|x|) is radial, then integrating in polar coordinates and taking into 
account that = (cr^- 1 ))- 1 / 2 we see that 

<^<o,i>= / F(x)<i, 1 (x)dx = c d f f(x)K /2 -\x)x d - 1 dx = c d (f,cf> d / 2 - 1 U> ld/2 _ 1 , 

JB d JO 

where Cd = (a(S d ^ 1 )) 1 ^ 2 . On the other hand, another integration in polar coordinates shows 
that for k > the Fourier coefficients (F, j, m ) vanish since then _L Y? = const. Thus the 

expansion of F in B d is in fact the expansion of its profile / with respect to the system {4>n 2 1 } 
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on (0,1). Similarly, the associated heat and subordinated semigroups are also related via the 
radial CELS6, clS stated below. 

Proposition 2.1. Let d > 1, v = d/2 — 1 and f G span{<^ : n > 1}. Then 

(0)o^x)=^' tt (/o P )(x), x £ B d , 

where f(x) = |x|. 

Proof. For / = (j)^ we have 

(TrX)o^x) = e- a -K(|x|) 

= Qe- tt M/-$5 x (x) = QT^^x) = 7^ a « o ^)(x), 

where q = (<t(S' c(_1 )) 1 / 2 . The conclusion follows. □ 

In fact the semigroups are related in the same way for more general functions /. This is 
confirmed by the relation between the corresponding integral kernels established below. 

Theorem 2.2. Let d > 1 and v = d/2-l. Then 

Gr(x,y)= [ g { t d) > a (x,yOdo-(Z), x,ye(0,l), t > 0, 
where x = (x, 0, . . . , 0) G B d . 

Proof. Let / G span{</>^}. Using Proposition 12.11 and then integrating in polar coordinates we 
get 

[ 1 Gr(\x\,y)f(y)y 2u+1 dy= [ g { t d) ' a (x,y)f(\y\)dy 

J0 JB d 

[ g { t d) ' a (x,yOda(C)f(y)y d - 1 dy. 

Since the subspace spanned by the <\> v n is dense in L 2 ((0, l),dfj, u ) we conclude the desired identity 
up to an exceptional set of y of measure 0. This set, however, must be empty because the kernels 
involved are continuous with respect to their arguments. □ 

Actually, each of the kernels G^' a (x,y) and ^^'"(x, y) is a jointly smooth function of t > 
and its arguments. This can be verified even directly, by term by term differentiation of the 
defining series, with the aid of basic bounds for d ntU , A n ^, J u (see [HI Section 2]) and the fact 
that |y^(£)l can be- dominated, uniformly in m and £, by a polynomial in k (see for instance [27} 
Chapter IV, Corollary 2.9 (b)]). We leave details to interested readers. 

3. Kernel estimates 

Heat kernels in various contexts were extensively investigated in the literature. In particular, 
it is well known that the heat kernel corresponding to the Dirichlet Laplacian in B d , d > 1, 
satisfies the bounds (cf. [17, (1.9.1)]) 

(2) <^ ),2 ( x 'y)^(i^ ex p(- ]£ ^ ! )' x 'y^ *>°- 

These bounds are perhaps most clear from the probabilistic point of view since C?^' 2 (x, y) is just 
the transition probability density of the time-scaled Brownian motion Bit killed upon leaving 
B d . However, finding sharp bounds for £/^' 2 (x,y) that describe precisely the interplay between 
x, y and t, and the boundary behavior, is a much more complicated matter. The complete and 
qualitatively sharp estimates are available only recently, for bounded C ' domains. The case of 
B d , obviously being geometrically simpler, does not seem to have been known earlier. 
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Theorem 3.1 ([T8J [TBI EQl [26] ). Let d > 1. Given T > 0, there exists a constant c > 1 such that 



(l-|x|)(l-|y| 



A 1 



1 



exp 



(3) 



<^ ),2 (x,y)< 



t d/2 

(i-M)(i 



c|x - y| 



A 1 



t d/2 



exp 



uniformly in x,y G i? d and < i < T. Moreover, ifT is chosen sufficiently large, then 
(4) ^ 2 (x,y)^(l-|x|)(l-|y|)exp(-tA^ /2 _ 1 ), x,y G B d , t>T. 

The upper bound in (J3j) as well as the large time behavior ^) has been known at least since 
the 1980s, cf. Davies and Simon |18| and Davies |16j . see also the related comments in |30| . 
The existence of the lower bound in ([3]) was shown by Zhang \30\ Theorem 1.1] under the 
assumption d > 3 and with an implicitly fixed T. Later these restrictions were removed by Song 
[261 Theorems 3.8 and 3.9]. Note that the upper estimate in (J3j> holds in fact for all t > 0, which 
is not the case of the lower bound. 

Sharp estimates for the subordinated kernels were found recently by Song [26] . 

Theorem 3.2 ((26j Theorem 4.7]). Let d > 1 and < a < 2. Given T > 0, we have 



(5) 



,(d),a 



(x,y) 



:i-N)(i-|y|) 

f2/a + | x _ y|2 



A 1 



(£2/a + | x _ y |2)(d+a)/2 ' 



uniformly in x, y G -B d and < i < T. 

Similarly as in ([3]), the upper bound in ([5]) holds in fact for all t > 0, which is not the case of 
the lower bound. The expected large time behavior for < a < 2 is of course 



,(d),a>. 



„ t f (x,y)^(l-|x|)(l-| y |)exp(-U^ /2 _ 1 ), x,y e B a , t>T, 

with T chosen sufficiently large. These bounds can be proved by a direct analysis of the series 
defining Q^' a (x, y), similar to that in the proof of Theorem 13 . 71 below. This, however, is beyond 
the scope of this paper. 

We now transfer the above bounds for Q^' a {pa, y) to the Fourier-Bessel setting on the interval 
(0, 1). More precisely, we shall prove the following. 



Theorem 3.3. Let v = d/2 — 1 for some d > 1. Given T > 0, there exists a constant c > 1 
such that 



( Xy y v - 1/2 (f a i 



v+l/2 



(l-a;)(l-y) 



<Gf(x,y)<(xy)-^ ^Al 



t 

xy 



A 1 



exp 



t 



v+l/2 



1 

{l-x)(l-y) 



yf 



A 1 



1 

Vt 



exp 



yf 



ct 



uniformly in x,y G (0, 1) and < f < T. Moreover, for sufficiently large T we have 

Gf(x,y)~(l- x )(l-y)exp(-t\l v ), x,y € (0,1), t>T. 
Theorem 3.4. Let a G (0, 2) and v = d/2 — 1 for some d > 1. Given T > 0, we have 



- {xy) 



-v-l/2 



(. 



xy 



A 1 



1^+1/2 



(l-x)(l-y) 

t 2/a + ( x _ y )2 



A 1 



t 



(t 2 /<* + (x-y) 2 )(<*+ 1 )/ 2 ' 



\t 2 / a + (x-y) 2 
uniformly in x,y G (0, 1) and < t < T. 

At the end of this section we will complement these results by deriving, directly and for all 
v > —1, sharp large time bounds for the kernels G^' a (x,y), < a < 2. 
To prove Theorem 13.31 we need the following. 
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Lemma 3.5. Let d > 1 and x,y,t,c> be fixed. Then 



exp 



where x = (x, 0, . . . , 0) € i/ = d/2 — 1 anc? ij, denotes the modified Bessel function of the first 
kind and order v. 

Proof. We first deal with the case d > 2. Observe that for £ £ £ 

|x-2/£| 2 \ ( x 2 + y 2 \ (Ixyi 



6XP ) = 6XP V cJ- ) 6XP 

Thus the integrand is a zonal function depending only on £i. In this case the integration over 
ga-l reduces to a one-dimensional integral against see for instance the proof of [271 Chapter 
IV, Corollary 2.16]. Precisely, we have 

J sd i ex P (^)dv(0 = v(S^) /_ 1 i ex P (^ei)(l-C 1 2 ) M/2 ^, 

where a{S d ~ 2 ) = 2^ d -^l 2 /T{[d - l)/2) is the surface area measure of S d ~ 2 . To evaluate the 
last integral we invoke Schlafli's Poisson type representation for the Bessel function I u (cf. [281 
Chapter III, Section 3-71, (9)]) 

Wz) = v^ + i/W ^* pM(1 - s2r ' /2 ' is - " > " 1/2 - 

Putting now all the facts together we arrive at the desired identity. 

When d = 1 the identity is verified directly, taking into account that (cf. |28} Chapter III, 
Section 3-71, (10)]) I_i/ 2 (*) = (2/(vrz)) 1 / 2 cosh z. □ 

Proof of Theorem \3.3[ The large time behavior is a consequence of Theorem 12.21 and Theorem 
13.11 To show the remaining estimates we use in addition Lemma [3. 51 This produces the bounds 



(xy) v 



(l-x)(l-y) 



< 



t 



A 1 



1 / x 2 + y 2 \ (2xy\ 



Gt (x,v) i$ (xy)' 



(l-x)(l-y) 



A 1 



r \ ct J \ ct / 



t 

Now the conclusion follows by applying the standard asymptotics for 7^, v > —1, 

Ij,^) ~ z", z — > + , /^(z) ~ z~ 1//2 exp(z), z — >■ oo, 

together with the fact that l v is continuous. □ 

In order to prove Theorem 13.41 we first analyze the behavior of a one- dimensional integral 
depending on several parameters. 

Lemma 3.6. Let 7 and n be such that 7 > r/ + 1 > 0. Then 
(l-s^ds 1 



< B < A < D. 



_! (D - Bs)(A - Bs)~' (D - B)A r i+ 1 (A - ' 
Proof. Observe that 

(l- s yds f 1 (l-s^ds f 1 (l-s)"ds 



Li (I? - - 5s)t J (D - Bs)(A - Bs)i J (D - Bs)(A - Bs)i ' 

so it is enough to analyze the last integral, which we further denote by X. It is convenient 
to consider the following three cases. Altogether, they cover all admissible configurations of 
A,B,D. 

Case 1: D > A > 2B. Since in this case A-Bs~A~A-B and D — Bs ~ D ~ D - B , for 
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s £ (0, 1), the conclusion is trivial. 

Case 2: D > 2B > A > B. Now D - Bs ~ D ~ D - B for s G (0, 1) and therefore 

1 ^(l-a^ds 



D - B J (A- Bs) 



Changing the variable u = 1 — s we get 



1 



it'' du 



D-B Jo (A-B + Bu)~< 

Splitting the last integral at the point u = (A — B)/B and estimating the resulting integrals 
separately we obtain 



(A-B)/B 



u r] du 



(A-B + Bu)t (A - B) 



(A-B)/B 



vP du 



Av+ l (A-B] 



■y-n-i 



and 



vP du 



< 



(A-B)/B 



(A-B + Bu)^ ~ Bf J (A _ B)/B 



u 



^ du < 



1 /A-B\v—r+i 



B7\ B 



An +1 (A- B)7-»7-i" 



This implies the desired conclusion. 

Case 3: 2B > D > A > B. Changing the variable u = 1 — s we get 

- 1 vP du 



o (D - B + Bu){A - B + Bup 

We now split the above integral at the points u\ = (A — B)/B and «2 = (D — B)/B and estimate 
the resulting integrals separately. We have 



(A-B)/B 



u v du 



1 



{D — B + Bu)(A — B + Bup (D-B)(A-Bp 

1 



(A-B)/B 



u n du 



Further, 



{D-B)/B 



u v du 



(D - B)A r >+ l (A - B)7-ij-i ' 

•(D-B)/B 



< 



1 



u 



(A _ B)/B (D-B + Bu)(A-B + Bup ~ (D - B)Brt J {A _ B)/B 

1 /A - By-v+i 
(D - B)B~< V B ) 
1 



^du 



< 



(D - B)A r >+ 1 {A- Bp-v-i 



Finally, 



vP du 



< 



1 



u 



{D-B) /B (D-B + Bu){A-B + Bu)~< ~ B^ J {D _ B)/B 

1 /D-B\v-~t 



'I-" 1 ' 1 du 



< 

< 



5 



1 



(D — B)A r > +l {A — 5)7-^-1 



The conclusion again follows. 



□ 



Proof of Theorem \3.4\ We assume d > 2; the case d = 1 is a simplified version of what follows. 
According to Theorem l2.2l we must integrate the right-hand side in ([5]) with x = (x, 0, . . . , 0) € B d 
and y replaced by y£, over S"^ -1 and with respect to da(£). To do that, we first use the relation 



A. NOWAK AND L. RONCAL 



z A 1 ~ z/(z + 1), z > 0, in order to switch to a comparable expression that is more suitable for 
the integration. This leads to 

f (l-x)(l-y) t 

l S d-i t 2 / a + |x ■ 



y£\ 2 + (1 - x)(l - y) {t 2 / a + |x - y ^)(d+a)/: 



Here |x — y£\ 2 = x 2 + y 2 — 2xy£i so, in fact, the integrated function is zonal. Thus (see the proof 
of Lemma [ 



Gt a {x,y) 



(l-x)(l-y)i 



/_i (t 2 / a + x 2 + y 2 + (l-x)(l-y)- 2xyii){t 2 / a + x 2 + y 2 - 2xy£,i) ( - d+a)/2 ' 

The last integral is handled by Lemma 13.61 applied with A = t 2 / a + x 2 + y 2 , B = 2xy, D = 
A + (1 — x)(l — y), 7 = (d + a)/2 and r/ = (d — 3)/2 = v — 1/2. Consequently, we obtain 

(l-x)(l-y)i 



(f/a + ( x _ y )2 + (1 _ a?)(l — J/))(t 2 /° + X 2 + y2) ^+1/2 (j2/ a + ( x _ y )2)(a+l)/2 ' 



This, combined with the relations x 2 + y 2 ~ (x — y) 2 + xy and z/(l + z) ~ z A 1, z > 0, finishes 
the proof. □ 

We close this section by describing the long time behavior of G^' a (x,y) for all u > — 1 and 
< a < 2. 

Theorem 3.7. Let v > — 1 and < a < 2. There exists T > such that 

G u t > a (x,y)~(l-x)(l-y)exp(-t\l„), x, ye (0,1), t >T. 
Moreover, the upper bound holds for an arbitrary fixed T > 0. 
Proof. We decompose, see ([I]), 

oo 
n=2 

To treat the first term here we note that 

<^(x)~(l-x), xG(0,l). 

Indeed, 4>\(x) is strictly positive for x G (0, 1), continuous for x G [0, 1] (see [15j (2.3)]), </>i(0) > 0, 
0i(l) = and Jjr^i(l) = — V2^i,u < 0. Here the last identity follows from the formula 



(6) 



7 't'n^.x) — x ^ A n i/ d n i/(A n i/x) ^ j/jy+i (A ra i/x) , 



which is a straightforward consequence of the differentiation rule for «/„, see |15|. (2.1)]. 
The proof will be finished once we show the bound 



£ 

n=2 



€(xWn(y) 



(1 - x)(l - y)e 



x,y€(0,l), t>T, 



for an arbitrary fixed T > and certain e > 0. To proceed, we first justify the bound 
(7) |C(x)| <(l-xK+ 2 , xG(0,l), n>l. 

By the Mean Value Theorem, for each x G (0, 1) there exists G (x, 1) such that 

d 



m X ) \ = m X ) -<(i)= (i-x) 
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To estimate the last derivative we combine © with the bounds given in |15t Section 2], see |15|. 
(1.3),(2.5),(2.6)], getting 

d " ' - Li n> h 9 e(0,l) 



<n(lV/ +1 / 2 )</ +2 



This gives ([7]). Now we may write 

oo 

^2e- tx ^\r n (xWn(y)\ < (l-x)(l-y)e- n 



oo 



-tc 



n=2 



n=2 



2i/+4 -i(A«„-A? -e) 



Choosing e such that < £ < Af — A" v and taking into account that A n>1/ ~ n, n — > oo (cf. 
[151 (2-6)]), we see that the last series is bounded uniformly in t > T, for any fixed T > 0. The 
conclusion follows. □ 

4. Comments and remarks 

Elementary cases. There are only two cases, v = ±1/2, when the Fourier-Bessel system has 
an explicit form, i.e. d n>v , J u and X n>u are all given explicitly. We have (see [21 Section 1]) 
K-i/2 = 7r(rt - 1/2), A„ il/2 = vrn, and 

(f>- 1/2 (x) = v^cos (vr(n - l/2)x), 



x 



l V2 



sin 7rnx . 



Consequently, the Poisson kernel G v t ' (x,y) can be computed when v 
basic trigonometric identities and the formulas (cf. [251 5.4.12 (1),(2)]) 



±1/2. Indeed, using 



£ e-^ sin(Bi) = - 



sini? 



we arrive at 



1 



2 cosh A — cos i? ' 



1 



cos(Bj) = - 



sinh j4 



2 cosh A — cos i? 



-1/2,1 



sinh ^ cos | (x — y) 



G 



1/2,1 



cosh 7rt 



1 



X ' y)= 2^ 



cos ir(x — y) 
sinh-7rt 



+ 



sinh y cos f ( x + y) 
cosh 7rt — cos 7r(x + y) ' 
sinh nt 



cosh 7rt — cos 7r(x — y) cosh 7r£ — cos 7r(x + y) 

Of course, sharp estimates for G^ 1//2,1 (x, y) can be obtained directly from these formulas, though 
it is not immediate. 

The heat kernel G^' 2 (x, y) is 'computable' as well for v = ±1/2, in the sense that the resulting 
expressions (in fact series) do not contain oscillations. For the simpler case v = 1/2 we have 



Gl /2 ' 2 (x,y) 



1 

xy 



E 



1 



4:7ft 



exp 



(x-y- 2jf 
At 



1 



'Ant 



exp 



{x + y- 2jf 
At 



This is essentially a well known Jacobi type identity. The series above represents the heat kernel 
in the setting of the system {V'nl with v = 1/2. Note that the corresponding differential operator 



is L v f(x) = x !/+1 / 2 L I ,((-)"^ 1/2 /)(x), and we have L 



1/2 



-A. The formula in question can 



be derived by solving the initial-value problem for the classical heat equation in the interval 
(0,1) with Dirichlet boundary conditions, see [H Chapter 3, Exercises 3.5-3.7]. Here the trick 
relies on extending a function / prescribing initial values on (0, 1) to an odd function on (—1, 1) 
and then considering its periodic extension f to R with period 2. The solution of the heat 
equation on R with initial values prescribed by /, which is given by convolving / with the 
Gauss- Weierstrass kernel, provides also the solution of our initial-value problem in the interval 
(0, 1). The convolution can be written in terms of an integral involving the above series, due to 
the symmetries of /. 

The case v = —1/2 is slightly more involved, but follows in the same spirit. Now L u = —A and 
one solves the initial- value problem for the heat equation in (0, 1) with the Neumann condition 
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at x = and the Dirichlet condition at x = 1 . The relevant extension / is obtained by extending 
/ to (—1,1) as an even function, then taking extension to (—1,3) that is antisymmetric with 
respect to x = 1 (i.e. satisfies f(x) = — f(2 — x) for x £ (—1,3)), and finally extending this 
function from (—1,3) to M as periodic with period 4. Solving the heat equation with initial 
values / and taking into account the symmetries of / one concludes that 

G^\x, y) = Y U= exp ( J*-V-V) 2 \ ^ _J_ ^ ( J* + V ~ W 





exp exp 



At 

(x + y- A(j + \)) 2 
At 



zbl /2 2 

The formulas for G t ' (x, y) are suitable to study short time behavior since for that only the 
terms corresponding to j = (the case of v = —1/2) or to j = 0, 1 (the case of v = 1/2) are 
essential. On the other hand, they indicate that the related analysis for general v > — 1 is far 

il /2 2 

from being trivial. We remark that the kernels G t ' (x, y) can be also expressed by Jacobi's 
elliptic theta functions 62,63 (see |25} p. 792] for the definitions), but this representation does 
not seem to be very useful for our purposes. 

Transference. Next we comment on the transference from the situation of Euclidean balls to the 
Fourier-Bessel setting on the interval (0,1). Since the semigroups in both settings are related 
via the mapping ip from Proposition 12. 1\ the same is true for many fundamental operators 
expressible through these semigroups. This concerns, for example, maximal operators of the 
semigroups, fractional integrals (potential operators) and Laplace transform type multipliers 
in both settings. Moreover, since Lebesgue measure in B d transforms to the measure dfj, v , 
v = d/2 — 1, when projecting via ip to (0, 1), mapping properties such as L p -boundedness and 
weak type estimates of operators defined in the context of Euclidean balls imply the same kind 
of mapping properties for the corresponding operators in the Fourier-Bessel setting on (0,1). 
Similar concepts of transference are well known also in other settings of classical orthogonal 
expansions, for instance in the direction Hermite — >Laguerre, see [21] and |23|. Section 5]. 

Rough upper bounds. Another issue to discuss are consequences of the rough estimate ([2]). 
Transferring it by means of Theorem 12.21 and Lemma 13.51 to the Fourier-Bessel setting on (0, 1) 
we get 

(8) 0<Gf (x,y) <W? +1/2 (x,y), x,y € (0,1), t > 0, 

where v = d/2 — 1, d > 1, and 

1 / x 2 + y 2 \ (xy 



W t \x,y) = (zv)-W*±e*p ( - 2-+!L)l A _ 1/a 



2/ 



is the heat kernel in the setting of continuous Fourier-Bessel expansions of type A > —1/2 on 
(0,oo) (the context of the Hankel transform), see [B] or [S]. The bounds ([S]) are easily explained 
from the probabilistic point of view. Indeed, W^(x,y) is the transition density for the time 
scaled Bessel process Xn on (0,oo), see Appendix I, Section 21] (here for A G (—1/2,1/2) 
the endpoint x = is assumed to be reflecting), and G t ' (x,y) is the transition density for 
X2t killed upon leaving (through x = 1) the interval (0, 1). Actually, this probabilistic argument 
works for all v > — 1 and justifies the following. 

Proposition 4.1. The bounds © hold for arbitrary v > — 1. 

This result has several interesting consequences. One of them is, see |24} Section 6], the fact 
that the semigroups {Tp a }t>o, v > — 1,0<q<2, are submarkovian symmetric diffusion 
semigroups (to get this for a < 2 we use the subordination principle). Since, in view of the prob- 
abilistic interpretation, both inequalities in ([8]) are strict, these semigroups are not Markovian. 
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Similar conclusions are valid for the corresponding semigroups in the setting of the system {V^} > 
with the restriction v E {—1/2} U [1/2, oo), see |24} Proposition 6.1]. 

Maximal operators. Further important consequences of Proposition 14.11 concern mapping 
properties of the maximal operators 

Tr' Q / = sup|r^ a /|. 

t>0 

Since T%' 2 is controlled by the analogous maximal operator W* +l ^ 2 in the Bessel setting on 
(0,oo) related to the measure d\i v , it inherits mapping properties of W* +l ^ 2 . The same is true 
in the context of the system {ip^} and the Bessel setting on (0, oo) related to Lebesgue measure. 
These observations allow to transmit positive parts of [61 Theorem 2.1] to the Fourier-Bessel 
setting. 

Theorem 4.2. Let v > — 1, 1 < p < oo, 5 £ R. Then the maximal operator T*' 2 , considered on 
the measure space ((0, 1), x & dx), has the following mapping properties: 

(a) T*' 2 is of strong type (p,p) if p > 1 and — 1 < 5 < {2v + 2)p — 1; 

(b) T*' 2 is of weak type (p,p) if — 1 < S < [2v + 2)p — 1 or 5 = 2v + 1; 

(c) T*' 2 is of restricted weak type (p,p) if —1 < 5 < [2v + 2)p — 1. 
Moreover, T*' 2 is of strong type (oo,oo). 

Analogous results, with appropriate adjustments, hold also in the setting of {ipn}, see [3 
Remark 3.2]. Boundedness properties with more general weights, and even in tensor product 
multi-dimensional Fourier-Bessel settings, can be deduced in a similar manner from the corre- 
sponding results in [5] and [9]; we leave details to interested readers. The above consequences 
complement and extend previous results on the heat semigroup maximal operator related to the 
system {V^} an d obtained in [20J and [12j . Note that, in view of the subordination principle, the 
abovementioned mapping properties are inherited by the maximal operators of the subordinated 
semigroups. Note also that Theorem 14.21 is not optimal in the sense that it does not take into 
account the boundary behavior at the right endpoint of (0, 1), which in comparison to the con- 
tinuous Bessel case is improved by the decay of the heat kernel, see Theorem 13. 31 Actually, this 
lack of optimality pertains to many other weighted results obtained so far for various operators 
in the Fourier-Bessel setting. 

The results of this paper shed some new light to the problems posed in [4] in which the heat 
semigroup comes into play. In particular, in connection with the question [U Q.6], it turns 
out that there is a pencil phenomenon associated with the heat semigroup maximal operator 
related to the system {V^}, and its nature is the same as in the case of continuous Fourier-Bessel 
expansions in the Lebesgue measure setting. Moreover, Theorem 13 . 3 1 enables to establish a sharp 
description of the pencil phenomenon in the Fourier-Bessel framework when the type index is 
half-integer. 

Final conjecture. We close the paper with the following natural conjecture. Proving (or 
disproving) it is an interesting and important open problem in the theory of Fourier-Bessel 
expansions. 

Conjecture 4.3. The estimates from Theorem \3.3\ and Theorem \3.4\ hold for all v > —1. 
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